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ABSTRACT 

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive 
method for observing the build-up of cosmological structure, which depends both on the ex- 
pansion rate of the Universe and our theory of gravity. In this paper we present a formalism 
for forecasting the constraints on the growth of structure which would arise in an idealized 
survey. This Fisher matrix based formalism can be used to study the power and aid in the 
design of future surveys. 
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1 INTRODUCTION 

The growth of large-scale structure, as revealed in the clustering of 
galaxies observed in large redshift surveys, has historically been 
one of our most important cosmological probes. This growth is 
driven by a competition between gravitational attraction and the 
expansion of space-time, allowing us to test our model of grav- 
ity and the expansion history of the Universe. Despite the fact 
that galaxy light doesn't faithfully trace the mass, even on large 
scales, galaxies are expected to act nearly as test particles within 
the cosmological matter flow. Thus the motions of galaxies carry 
an imprint of the rate of growth of large-scale structure and al- 
lows us to probe both dark energy and test General Relativity (e.g . 
. [ Jain & Zhanj2008l:ISong & Kovamj|2008l:ISong & Percivall2008l ; 
IPerciv al & Whit el2008l ; lMcDonald & Selialj2008t for recent sttid- 
ies). 

This measurement of the growth of struc ture relies on r edshift- 
space distortions seen in galaxy surveys jKaised [l987h . Even 
though we expect the clustering of galaxies in real space to have 
no preferred direction, galaxy maps produced by estimating dis- 
tances from redshifts obtained in spectroscopic surveys reveal 
an anisotropic galaxy distribution. The anisotropics arise because 
galaxy recession velocities, from which distances are inferred, in- 
clude components from both the Hubble flo w and peculiar veloc- 
ities driven by the clustering of matter (see lHam ilton 1998j, for a 
review). Measurements of the anisotropies allow constraints to be 
placed on the rate of growth of clustering. 

Ever larger surveys have provided ever tighter constraints. 
Anal yses using the 2-d egree Field Galaxy Redshift Survey (2dF- 
GRS; IColless et al.l2003h have measured redshift-sp ace distortions 
in both the correlation function ( Peacock et al. 2001 ; Hawkins et al. 
2003) and power spectrum ( Percival et a lj 120041). Using the 
Sloan Digital Sky Survey (SPSS: lYorket alJ l200d) , redshift- 
spac e distortions have also been measured in the correlation func- 
tion ( Zehavi et al ] |2005l : lOkumura et aT]|2008l ; ICabre & Gazta naga 
2008), and using an Eigenmode decomposition to separate real 



and redshift-space effects dTegmark et ai]|2004l2006|). These stud 



ies were recently extended to z — 1 I Guzzo et alj 2007 ) us- 
ing the VIMOS-V LT Deep Survey (VVDS: iLe Fevre et alj|200l 
Garilli et alj|2008h . In addition to measuring clustering growth at 
z = 0.8, this work has emphasized the importance of using large- 
scale peculiar velocities for constraining models of cosmic ac- 
celeration. Current constraints on the growth rate are at the sev- 
eral tens of percent le vel (e.g. iNesseris & Perivo laropoulosj |2008l : 
ISong & Per cival 2008), but observational progress is rapid. 

In the next section we shall outline the formalism for forecast- 
ing constraints on cosmological quantities from measurements of 
redshift space distortions, and compare it with previous forecasts. 
We begin with the simplest model and then investigate various re- 
finements. We finish in ^4]with a discussion of future directions. 
For illustration we shall assume a fiducial ACDM cosmology with 
Q m = 0.25, h = 0.72, n = 0.97 and <x 8 = 0.8 (in good agreement 
with a variety of observations) when computing specific predictions 
for future surveys. 



2 THE FISHER MATRIX 

The Fisher matrix provides a method for determining the sensitivity 
of a particular measurement to a set of parameters and has been ex- 
tensively used in cosmological forecasting and optimization. Here 
we adapt this methodology to our particular problem. 



2.1 The simplest case 

Under the assumption that the density field has Gaussian statistics 
and uncorrected Fo urier modes, the Fish er matrix for a set of pa- 
rameters [pi] is (e.g. iTegmarketalJ 19981) 
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where P is the power spectrum and the mode cou nting is deter- 
mined by the effective volume dFeldmanetalJ 19941) 



(2) 



which depends on the geometric volume of the survey, Vo, and the 
number density, h, of the tracer. If h is high enough then V e ff - Vo. 
The constraints are dominated by regions where hP > 1, so it is safe 
to neglect the higher order (in fT 1 ) terms which arise assuming that 

alaxies are a Poisson s ample of the underlying density fluctuations 

Meiksin& Whitdl 1999h . 

The simplest model for the observed galaxy distribution 
is a linear, deterministic, and scale-independent galax y bias, 
with redshift space distortions due to super-cluster infall iKaiseil 
and no observational non-idealities. In this case /^ b s oc 
(b + //J 2 ) P\m(k) where Py m is the linear theory mass power spec- 
trum in real space, b is the bias and fj the angle to the line-of-sight. 
The quantity of most interest here is / = d In D/d In a, the logarith- 
mic derivative of the linear growth rate, D(z), with respect to the 
scale factor a = (1 + z) -1 . In general relativity / « Q. mat (zf 6 (e.g. 
|Peebleslll98Ch . while in modifie d gravity mode ls it can be smaller 
by tens of percent (e.g. lSong & Percivalll2008t figure 1). Redshift 
space distortions allow us to constrain / times the normalization of 
the power spectrum (e.g. f(z)cr i (z)), or dD/dlna. The derivatives 
in Eq. (TJ are particularly simple 
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independent of the shape of the linear theory power spectrum, and 
hence of the spectral index and transfer function. Since we hold the 
normalization of the power spectrum fixed for these derivatives, the 
fractional error on f(z)crg(z) is equal to that on / in our formalism. 
The errors on b and / depend sensitively on the maximum k in the 
integral of Eq. |QJ. Since we are using linear theory we choose to 
cut the integral off at k = 0. 1 h Mpc 1 for our fiducial cosm ology. 
This is close to the value at which IPerciva l & White ( 200§) saw 
departures from linear theory. 

The bias and / turn out to be anti-correlated, with a correla- 
tion coefficient of 70 - 75%, depending on the precise sample. We 
marginalize over b by first inverting the Fisher matrix to find the co- 

/ ,\l/2 

variance matrix and hence the error /, specifically Sf = \F j . 
Hereafter we shall drop the parentheses and write for the // 
component of F -1 . 

Forecasts on Sf If in this model can be regarded as an "op- 
timistic" estimate of the reach of future observations, and the re- 
sults are quite encouraging. For a 10(/r'Gpc) 3 survejQ at z = 
an unbiased population of tracers with hP » 1 would achieve 
Sflf = 1.6%, with the error scaling as Vq 1/2 . The constraint 
weakens with increasing bias, being 3% for b = 2 and 6% for 
b = 4. The constraint also weakens as we reduce h. For b = 1 
and n = 10~ 3 h 3 Mpc~ 3 the constraint is 1.8%, increasing to 1.9% 
for h = 4 x 1CT 4 h 3 Mpc~ 3 , 2% for n = 2 x 1(T 4 h 3 Mpc~ 3 and 3% 
for n = 10- 4 A 3 Mpc" 3 . The shot noise, nP(0.1 h Mpc" 1 ) =* 1 for 
h = 2 x 10~ 4 h? Mpc~ 3 , which explains the rapid increase in Sf for 
h larger than this. 

Conversely, increasing & max to 0.2/iMpc" 1 reduces the error 
on the b = 1, hP » 1 case to 0.6%. Moving to higher red- 
shift, keeping the bias fixed, makes the constraint stronger as f/b 



1 Out to z = 1 over 10,000 square degrees the (comoving) volume in our 
fiducial cosmology is 14(/i~'Gpc) 3 . 



is increased. By z = 1, / has increased to 0.83 from 0.44 and 
Sflf - 1% for our fiducial b = 1, hP » 1 example. Note, however, 
that at higher z the effects of shot noise would typically be larger. 

These cons traints can be compared to the forecasts in 
Guzzo et al. I d2007h who present a fitting function for the relative 
error on/3 = f/b of 



(4) 
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where h is measured in h 3 Mpc~ 3 and V in h~ 3 Mpc 3 . Note that both 
forecasts agree on the scaling with volume, but the above scales ap- 
proximately as the inverse square root of the total number of galax- 
ies in the survey and is independent of the bias. Taking into account 
the correlation between the constraints on b and / our forecast con- 
straint is 

S j=b- l [p>F- b l-2pF- b } + F-/ f ] m (5) 

Comparing our forecasts to this scaling we find rel atively good 
agreem ent for b = 1 and h - 10~ 4 h 3 Mpc -3 , but the lGuzzo et al.l 
( 2007) scaling predicts much better constraints for higher number 
density or more biased samples. 



2.2 Beyond linear theory 

Of course we do not expect the simple linear theory result with 
super-cluster infall to be a perfect description of redshift space dis- 
tortions on all scales. Comparison with N-body simulations sug- 
gests that halos do closely follow the matter velocity field and 
the major modification to the simple model at low k is in the 
quadrupole, with an additional effect coming from the generation 
of multipoles higher than 4. By introducing more freedom into the 
model we will increase our ability to describe the extra physics act- 
ing, an d simultaneously b e gin to degrade our sensitivity to /. 

In IPercival & White! d2008h it was shown that a streaming 
model with a Gaussian small-scale velocity provided an adequate 
fit to N-body simulations to k m 0. l/iMpc~'. Under these as- 
sumptions one model for the redshift space, galaxy power spectrum 
could be 

P Bbs (k, n) = {b + f//f P (k)e- k2 ^ 2 (6) 

where P represents the mass power spectrum in real space and <r- 
is to be regarded as a fit parameter which encompasses a variety 
of violations of the traditional analysis. We can additionally model 
inaccuracies in the observed redshifts by a line-of-sight smearing 
of the structure. In the limit that this smearing is Gaussian it can be 
absorbed into o%. 

In this situation, the logarithmic derivatives with respect to 
b and / are unchanged and the new derivative required is simply 
3 In Pi 'da 2 . = -k 2 fi 2 . We now marginalize over <x, in addition to b 
before reporting the constraints on /. 

Many of the trends with b and h in our simple model also hold 
for this extended model. For our fiducial 10(/r'Gpc) 3 volume the 
constraint from our z = 0, unbiased tracers with hP » 1 weakens 
from Sflf - 1.6% to 3.2% when marginalizing over cr z . The fidu- 
cial value of <X; has little impact for nP » 1. However when hP = 1 
the error is increased from 3% to 4% to 20% as cr z is increased from 
to 10 /r'Mpc to 100 h~ [ Mpc for a sample with b = 1. 

Another alternative is to model the small-scale suppression 
with a Lorentzian, which provides a better fit at higher k and is a 
good match to the superposition of Gaussians of different widths 
from halos of different masses JWhitdl200lh . The two agree to 
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Figure 1. The correlation coefficient, r(k), between the density and velocity 
divergence of the dark matter in an N-body simulation of a ACDM cosmol- 
ogy, with the same cosmological parameters as our fiducial model. 




lowest order in kcr z , and which is the regime of most interest 
here, with the Gaussian matching the results of N -body simulations 
at sm all k slightly better than the exponential jPercival & White! 
2008). Changing the form from Gaussian to exponential makes a 
negligible change in our forecasts. 



2.3 Mode by mode 

The forecasts above all made quite strong assumptions about the 
relationship between the velocity and density power spectra, as- 
sumptions which are only known to be true for quasi-linear scales 
within the context of General Relativity. The parameters j in our 
Fisher matrix (Eq. QJ don't have to be cosmological parameters 
however. We can fit directly for the three independent power spectra 
(the density-density, velocity-velocity and density- velocity spectra) 
rather than assuming that they are related by a specific functional 
form (e.g. iTegmark et al]|2004r) . Such constraints would be appli- 
cable to a wide range of theories including e.g., interacting dark 
energy, clustered dark energy or f(R) gravity. 



Figure 2. (Top) Pm(ki,Zj) in redshift bins of width Az = 0.2 at Zj = 0.5 
(solid) and 1.5 (dashed) from a half-sky survey with b = 1,5 and n = 5 X 
10~ 3 h 3 Mpc~ 3 . (Bottom) The fractional error on P@@(kj,Zj) in the same 
bins weighting modes with o",i, = 0.1 (see text). 



the two spectra themselves, in bins of k and z: 



1 



P obs (ki,iu,Zj) 



din P obs (k h n,Zj) 
dP gg (ki,Zj) 

din P obs ,(ki,Zj) 
dP && (ki,Zj) P ohs (ki,n,Zj) 

and the variance of P @@ (k h Zj) is given by 

<T[Pee(k„Zj)] = F^(k„ Zj ) 
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where F^ 2 ' (kt, Zj) is 22-component of the inverse matrix of F a p. The 
constraint on P &0 (kj, Zj) in 2 redshift bins each of Az = 0.2 is plot- 
ted in Fig. [2] for a half-sky survey with h = 5 x 10~ 3 h 3 Mpc~ 3 and 
b = 1.5. 



2.3.1 Correlation of 6 and ® 

N-body simulations of ACDM cosmologies show that the density 
and velocity divergence are highly correlated for k < 0.1 /i Mpc~' 
(see Fig.[]} so we will begin by making the assumption that the den- 
sity and velocities are perfectly correlated (to be relaxed in £|2.3-2t . 
Then the density-velocity cross-spectrum becomes the geometric 
mean of the two auto-spectra and we have only two free functions. 
If we write for the velocity divergence in units of aH the power 
spectrum becomes 



Pob S (k,n,z) = P gg (k,z) + 2n 2 [p gg (k,z)P 00 (k,z)] 



1/2 



+ i/P &e (k,z) 



(7) 



where P gg denotes the usual galaxy density aut o-spectrum and we 
have assumed that small-scale ("finger of god"; Ijacksonll 19721) ef- 
fects have been cleanly removed by e.g. finger of god compression. 
The parameters in the Fisher matrix, Eq. l[T}, are now the values of 



2.3.2 Decorrelation of 6 and © 

The assumption of tight correlation between 5 and is a reasonable 
one (densities grow where flows converge and velocities are high 
where mass concentrations cause a large gravitational potential) but 
is not required. We can extend the formalism above by allowing the 
cross-correlation coefficient, 

r(k) = 80 , (10) 

^P && {k)P gg (k) 

to differ from unity. The power spectrum can now be written in 
terms of 3 free functions (P gg , Pqq, r) as 

P b, = (P gg + 2ir r(k) ^P gg P 00 + fi 4 P @@ ) G FoG (k,^; cr z ) (11) 

where Gf g is a Gaussian describing the decrease in power due to 
virial motions and the derivatives are given by 
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(12) 



We find that allowing r(k) to be completely free degrades 
the constraint on P@@, and eventually /, significantly, until it is 
equivalent to simply measuring the fi 4 component in Eq. Jilt . 
To strengthen the constraint requires prior information about r(k), 
which can in principle be obtained from simulations or perturba- 
tion theory calculations of structure formation in modified gravity 
models. As an illustrative example, if we assume a prior 



, r (r) = 1 - r(k) 



(13) 



where the error on r is equal to its deviation from unity (using 
the fiducial model with r measured from N-body simulations as 
in Fig. [TJ we find the constraint on P @@ is almost the same as we 
obtained before. 

Also note that in this analysis we can mitigate our uncertainty 
in the form of the small-scale redshift space distortion by_ down- 
weighting modes for which the fingers-of-god (JacksoH I 19721) are 
expected to be large. The residual uncertainty after the weighting is 



1 r V d 3 k 

2 J (2/r) 3 
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dpi )\ dPj 



V eff (k,n) w FoG (k,fi) 



(14) 



where the weight function wp oG (k,fj.) could, for example, be given 

by 



w FoG (fc,yu) = exp 



(G F , 



1) 



(15) 



where Gf g is the finger-of-god suppression factor and tr th is a 
threshold value indicating our confidence in the accuracy of the 
FoG model. 



2.4 Multiple populations 

Until now we have implicitly assumed that we are dealing with a 
single population of objects. However galaxies come in a variety of 
sizes, luminosities, masses and types which exhibit different clus- 
tering patterns but all of which are expected to r espon d to the same 
large-scale velocity field. McDonald & Seljak (2008) pointed out 
recently that this allows, in principle, for significant gains in deter- 
mination of the growth of structure. In fact, in the limit of Gaussian 
statistics, perfectly deterministic bias and infinitely dense tracers, 
one can measure the velocity power spectrum limited only by the 
total number of modes in the survey in all directions. 

To include multiple populations in the Fis her matrix approach, 
there are two obvious ways of proceeding. McDonald & Seljak 
d2008h assumed that the densities, <5, for 1 < i < N, are the mea- 
sured quantities and built a covariance matrix in terms of the power 
spectra, (<5,(5;), where a superscript denotes a measured quantity 
that includes a noise term. An alternative and complementary ap- 
proach is to extend the analysis presented in §2. ll assuming that the 
power spectra are the measured quantities. For Gaussian fluctua- 
tions, in which all of the cosmological information is encoded in 
the power spectrum, these approaches turn out to be equivalent 
We develop this second approach here. 
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Figure 3. The fractional error on fiz)o'g(z) arising from a 10(/r'Gpc) 3 
survey at z. = populated with two types of galaxies. The first population is 
held fixed with b[ = 1 and hi = 10~ 2 h 3 Mpc~ 2 , i.e. nP » 1. The second 
population has b2 = 1.4 (solid), bi = 2 (dashed) or \>2 = 4 (dotted) and the 
constraint is plotted vs. «2. All else being equal, the fractional constraints 
would be tighter at higher z where / is larger. 



2.4. 1 The Fisher matrix 

The Fisher matrix for this problem is a simple generalization of 
Eq. l[T) but now there are N(N + l)/2 power spectra for A' galaxy 
populations. For example, in the simplest case of two populations 
there are 3 measured power spectra, which on large scales are esti- 
mates of 



(16) 



where i and j run over 1 and 2 and P\o = P%\ ■ 

To calculate the Fisher matrix for multiple samples we need to 
sum over [N(N + l)/2] 2 elements of the inverse covariance matrix 



z 

XI' 



rvofk (dP x \ j tdPy 

J (2n) 3 [dp~r xr [dp-j 



(17) 



where we denote a pair of galaxy indices by X or Y. Note that this 
reduces to Eq. Q} for a single population. In order to calculate the 
Fisher matrix, we need to determine the covariance matrix and the 
derivatives of the power spectra with respect to the cosmological 
parameters of choice. 



2.4.2 Calculating the covariance matrix 

If we assume that the bias is deterministic and that the shot- 
noise can be treated as an (uncorrelated) Gaussian noise the co- 
variance matrix for the power spectra is straightforward to com- 
pute. Similar results for the covariance matrix of quadratic com- 
binations of Gaussian fields have been determined previously 
when c onsidering CMB temperature and polarization power sp ec- 
tra (e.g. lZaldarriaga & SeliaMll997l ; lKamionkowski et al.lll997l) or 
the problem of combin ing density and velocity power spectra (e.g. 
iBurkev & Tavlorll2004l) . Our problem is slightly more general, in 
that we need to consider additional combinations of power spectra, 
but similar in spirit. 

If we define N a = [1 + l/(hP aa )], the diagonal terms in the 
covariance matrix are 



We verified this by explicit numerical computation of the Fisher matrices. 



(18) 



Redshift-space distortions 



(Cabab) ~ P a b + P aoP bbN a N b , 

where a 4- b, and the off-diagonal terms are calculated from 
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where a 4 b 4- c 4- d. These formulae are complete with the rela- 
tions P a t = Pta and (Cxy) = (Cyx)- The off-diagonal terms given 
in Eq. &221 do not occur in the CMB example as there is only a 
single non-zero cross power there. However, all covariance matrix 
elements can be calculated using the same standard procedure (see 
Appendix lAl). 



2.4.3 Calculating the derivatives 

If the parameters that we wish to constrain are b a and / then 
^ = [(b b + f^ + (b a +f^]P, m , 

d -j± = [(b b +f^ 2 )+(b a +f^)yp lm 

where 6"£ and 6 bc are Kronecker Ss. This completes the input that 
we need for the Fisher matrix, Eq. j!7\ . 

We can include a parametrized line-of-sight smearing by mul- 
tiplying the power spectra by e.g. exp[-(l/2)A: 2 /j 2 (cr^ + cr 2 b )]. The 
derivatives are multiplied by the same factor and there is an addi- 
tional set 



(24) 
(25) 
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dPqb _ 1 ; ,2 , ,2 / xac . x bc\ n 

— = --kfj [6 K +8 K )P a , 



For the mode-by-mode parametrization developed in §23] we 
can use the logarithmic derivatives in Eq. d 12b . multiplied by 



2.4.4 Results 

We confirm the finding of McDonald & Seliakl d2008l) that using 
multiple populations can result in significant gains in constraining 
power. For example Fig.|3]shows the fractional error on / decreases 
by a factor of 2 - 3 if a second population is simultaneously used to 
provide constraints. At fixed number density the gain is higher the 
more biased is the second sample, and the constraint is weakened 
as the bias of the first sample is increased. Thus we would like 
to find a two samples with very different clustering properties but 
reasonable number densities. 

As long as the line-of-sight dispersion, <r,, is not large the 
marginalization has little effect on the total error. One does, how- 
ever, prefer slightly higher hP when marginalizing over tr, than 
when keeping it fixed. 

The gains saturate quickly when using more than two samples. 
In fact if the total number of objects observed is to be held fixed, it 
is better to increase the number densities of the lowest and highest 
biased sets rather than include an intermediately biased sample at 
the expense of lower number densities for all samples. 

Within the deterministic bias model, splitting into multiple 
populations does not affect constraints on the overall large-scale 
power spectrum shape: here we are always limited by the total num- 
ber of modes in the sample. A bias model can be used to weight 
galaxies of different bias, allowing for their different clustering 



BOSS 

WFMOS (1) 
WFMOS (2) 



3 0.1<z<0.7 1.5 
5 0.5<z<1.3 2.0 
5 2.3 < z < 3.3 0.6 



EUCLID/JDEM 50 0.1 < z < 2.0 500 



Table 1. Fiducial parameters adopted as indicative of various planned or 
ongoing surveys. N„z\ is given in units of 10 6 , e.g. BOSS has 1.5 million 
galaxies, and h in units of 10~ 4 ft 3 Mpc~ 3 . For the survey with the proposed 
WFMOS instrument, we assume tha t this is split into low (1) and high (2) 
redshift components as proposed in Glazebrook et al. 12005). We assume 
that each survey covers a fixed fraction of the sky, so the volume within any 
redshift interval is completely determined by these parameters. 
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Figure 4. The fractional error on /(z)o"8 (z) in bins of width Az = 0. 1 , arising 
from fiducial surveys with parameters given in Table fJJ We consider all 
galaxies in a single bin (solid lines), and split into 4 bins as a function of 
bias (d ashed lines). The existing constraints, as collected in lSong & Percivaj 
1 2008) and with the addition of da Angela et al. ( 2008), are shown as solid 
squares (see text). 



strengths, in order to optimally calculate the overall power spec- 
trum shape jPercival et alj|2004h . Any cosmological benefit from 
splitting into multiple samples will therefore arise through better 
constraints on f(z)crg(z). 



3 PREDICTIONS FOR FUTURE SURVEYS 

In this section we apply our Fisher matrix formalism to 3 concepts 
for future spectroscopic survey^ with fiducial parameters given in 
Table[T] We assume a tight prior on small-scale velocity dispersion. 
The galaxy bias is one of the hardest parameters to predict for fu- 
ture surveys, so we have adopted a conservative approach here. We 
assume that redshift zero galaxy bias is sampled from a uniform 
distribution with 1 < b < 2. The bias evolves with redshift such 
that the galaxy clustering amplitude is constant. For all surveys, we 
assume that we can use all modes with k < 0.075 h Mp c~' at z = 0, 
and that this limit evolves with redshift according to the Smith et al. 
d2003h prescription for k n i. This assumption is deserving of further 
investigation in N-body simulations. 



3 The HETDEX experiment fall et aljpiooi) has constraints similar to 
WFMOS(2) and we have not plotted it to avoid crowding. 
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BOSS (z=0.35) 
WFMOS(l) (z=0.65) 
WFM0S(2) (z=1.65) 
EUCLID (z=l) 




Figure 5. The fractional error on Pqq using the experimental parameters 
in Table [T] We show only a single, representative bin in redshift for each 
experiment, to avoid clutter. 



Predictions for the error on /(z)<x 8 (z) are shown in Fig. [4] Re- 
sults are presented either assuming that all galaxies are analyzed in 
a single bin, or are split according to galaxy bias. For surveys with a 
large number density, such as that proposed for the EUCLID/JDEM 
concept, splitting into bins with different galaxy bias can signifi- 
cantly reduce the expected errors, as we saw in i ]2.4| The galaxy 
sampling density proposed for the BOSS and WFMOS surveys is 
lower, and we gain less from the multiple-sample approach. 

It is useful to ask how these forecasts depend on the input as- 
sumptions. As an illustration, if we decrease the number density by 
a factor of 1.5 the BOSS, WFMOS(l) and unsplit Euclid/JDEM re- 
sults are largely unchanged. The error on WFMOS(2) increases by 
30% while the split Euclid/JDEM limit increases by 10%. Decreas- 
ing the maximum k from 0.075 h Mpc~' to 0.05 hMpc~' at z = 
(with the same scaling to higher z) all of the constraints become 
weaker. The limit from BOSS increases by ~ 40%, WFMOS(l) by 
~ 75%, WFMOS(2) by ~ 15% and Euclid/JDEM by ~ 85% for the 
unsplit case and 60% for the split case. 

We also i nclude current constraints, as collected by 
ISong & PercTvaH J2008h and with the addtion of 2SLAQ, in Fig.|] 
These are 6(f(T & )/(f cr s ) = 0.12 at z = 0.12 from the 2dFGRS 
jPercivaletal]|2004h , S(f( T S )/(fcr & ) = 0.11 a t z = 0.35 from 



the SDSS LRG catalogue iTegmark et al 

0.37 at z = 0.77 from the VVDS JGuzzo et al.1 |2007|) and 



2006), S(f<T*)/(fc r x ) = 



1.5 from 2SLAO jda Angela et al.1 



6(fc£)/(f<r t ) = 0.21 at z 
2008). As can be seen, the next generation of spectroscopic galaxy 
surveys will provide an order of magnitude increase in the available 
cosmological constraints from redshift-space distortions at z > 0.1. 

The direct constraint on Pq@ for our futuristic surveys is 
shown in Fig. [5] for some representative bins in redshift. We expect 
to be able to place tight constraints on the velocity power spectrum 
near k ^ OA hMpc~' with future experiments. 



4 CONCLUSIONS 

Observations of redshift-space distortions in spectroscopic galaxy 
surveys offer a powerful way to measure the large-scale velocity 
field, which in turn provides a sensitive test of both the expan- 
sion rate of the Universe and our theory of gravity. We have de- 
veloped a Fisher matrix formalism which allows forecasting of the 
constraints future, idealized, surveys would be able to place on the 
linear growth rate, /(z)cr 8 (z) oc dD/d\na, and shown that they are 
potentially highl y constraini ng, though not as constraining as the 
scaling of lGuzzo et alj <2007h predicts. 

We have developed the Fisher matrix exposition in multi- 
ple levels of sophistication and realism, assuming strict functional 
forms for the power spectra or allowing them to float freely. As ex- 
pected the constraints are tightest when theoretical investigations 
can provide good priors for the form and range of parameters, but 
even relatively conservative assumptions suggest that percent level 
measurements of / should be possible with future surveys. Fur- 
ther work on understanding the correlation of velocity and density 
fields, scale-dependent bias and non-linear effects could pay big 
dividends. 

As pointed out by iMcDonald & Seliakl fe008h . using multi- 
ple populations of galaxies can tighten the constraint on /(z)<xg(z) 
(though it does not improve measurement of the shape of P@@). 
We show that this can be naturally incorporated into our formal- 
ism. The largest improvement comes when disjoint samples with 
a large difference in bias, both having a high number density 
[hP(k O.lfcMpc -1 ) » 1] are used. The ultimate limit to this 
method will come from stochasticity in the biasing of galaxies, and 
which types of galaxies minimize this effect on which scales is an 
important avenue for further investigation. 

We have made a number of simplifications in this analysis 
which it will be important to address in future work. First we have 
assumed that the large-scale velocity field of the galaxies is that of 
the matter. Ultimately our ability to model any velocity bias will set 
a lower limit on what can be achieved. It is important to note that we 
are limited by how accurately the halo centers follow the mass ve- 
locity field "on large scales", rather than a bias within the halos. For 
the forme r, simulations suggest that halos do tend t o trace the mass 
very well fauff et al.l2007h|Percival & Whitekooll) . Current obser- 
vational constraints on the latter from modeling c lusters are consis - 
tent with no velocity bias at the 10% level (e.g. Sodr e et al.1 1989). 
Simulations suggest little or no velocity bias for the majority of 
"galaxies" in dark matter ( Springel et al. 2001; Faltenbacher et al. 
2006) and hydrodynamic (Berlind et al. 2003) simulations at the 
same level. Investigations of these phenomena in simulations of the 
standard cosmology and with alternative theories of gravity will be 
very valuable. 

A code to compute the Fisher ma- 
trix given survey parameters is available at 

|http : //mwhite . berkeley . edu/Redshif t| 
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APPENDIX A: THE OFF-DIAGONAL POWER 
SPECTRUM COMPONENTS 

In this section, we derive formulae for {C aa bc}, and (Cabac), as given 
in Eqns.l21l& l22l Similar derivations fo r CMB power spectra were 
presented in lKamionkowski et alj ( 119971) . 

Suppose that, in a particular experiment, we have M indepen- 
dent complex samples S"\ with 1 < m < M, drawn from multivari- 
ate Gaussian distribution with P ab = l/M ZmW^i )• O ur estimate 
of the covariance between power spectrum measurements from this 
experiment is 

(PabPcd) = ^ 2<«<5"'X > (AD 

mjn' 

To proceed, we split this sum into terms with m = m' and m ± 
m'. Where m = m', we can use the standard result for the 4-order 
moments of multivariate Gaussian random variables that, if x a are 
real and Gaussian distributed, the expectation 

E[x a x b x c x d ] = (x a x h ){x c x d ) + {x a x d ){x b x c ) + (x a x c )(x b x d ). (A2) 

For the component of Eq. lAll where mj=m', we can easily decom- 
pose into 2-order moments. For {P aa P b c) this procedure gives 

(PM = ^£[<IOXT><5™"' + 

m,m' 

{\6 m M''S'"'){l-5T')\ (A3) 

= — \ PaaPbc + 2P ab P ac \ + P a „P bc - —PaaPbc (A4) 
= — PabPac + PaaPbc, (A5) 

M 

so (Caabc) = 2P a bPac divided by the number of modes. 

To calculate C a b ac , note that this procedure also gives that 

(PabPac) = l/M[P a bPac + PaaPbc] + PabPac, SO (C a bac) = PaaPbc + 

PabPac, divided by the number of modes. The other terms in the 
covariance matrix can be calculated using the same methodology. 



